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TSALLIS STATISTICS AS A TOOL FOR STUDYING INTERSTELLAR TURBULENCE 
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ABSTRACT 

We used magnetohydrodynamic (MHD) simulations of interstellar turbulence to study the probability distri- 
bution functions (PDFs) of increments of density, velocity, and magnetic field. We found that the PDFs are well 
described by a Tsallis distribution, following the same general trends found in solar wind and Electron MHD 
studies. We found that the PDFs of density are very different in subsonic and supersonic turbulence. In order to 
extend this work to ISM observations we studied maps of column density obtained from 3D MHD simulations. 
From the column density maps we found the parameters that fit to Tsallis distributions and demonstrated that 
these parameters vary with the sonic and Alfven Mach numbers of turbulence. This opens avenues for using 
Tsallis distributions to study the dynamical and perhaps magnetic states of interstellar gas. 
Subject headings: ISM: general — MHD — turbulence 



1. INTRODUCTION 

Turbulence is known to be a crucial ingredient in many 
astrophysical phenomena that take place in the interstellar 
medium (ISM): such as star formation, cosmic ray dispersion, 
formation and evolution of the (also very important) mag- 
netic field, and virtually every transport process (see for in- 



stance the rev i ews of Elmegre en & Scal o 2004 ; Mac Low & 
Klessen|2004"l|Ballesteros-Paredes et al.|2007HMcKee & Os- 
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triker||2607 and references therein). This has brought many 
researchers to the study of turbulence, from both observa- 
tional and theoretical perspectives. 

From the observational point of view there are several tech- 
niques aimed to study propertie s of ISM turbulence. Scintil- 
lation studies (see for in stance |Narayan & Goodman||1989[ 
Span gler & GwinnJT 990) have been very successful to char 



acterize ISM turbulence, however they are limited to den- 
sity fluctuations in only ionized media. For neutral media 
the density fluctuati ons can also be obtained from column 
density maps (see Padoan et al. 2003). Other avenue of re- 
search is provided by radio spectroscopic observations , using 
for instance the centroids of spectral line s ( von Hoerner| 1951] 
Munch|19581|Dickman & Kleiner) 1 985] |Kleiner & Dickman l 
1985||Q'Dell & Castaneda||1987[fFalgarone et al.|1994[|Mi-| 
esch & Bally|1994| |Miesch & Scalo 199 5||Lis et al.|1998j|MTi 
esch et al. 1 1 999| ), or linewidth maps ( |Larson| 1 98 1 1 1 1 992| |Scalo| 
1984||1987[ ). Spectroscopic observations have the advantage 



of containing information about the underlying turbulent ve- 
locity field. However, they are also sensitive to fluctuations in 
density, and the separation of the two contributions alone has 
proved to be a difficult problem (see Lazarian 2006b). 

Much of the effort to relate observations with models of 
turbulence has focused on obtaining the spectral index (log- 
log slope of the power-spectrum) of density and/or velocity. 
To this purpose several methods have been developed and 
tested through numerical simulations. Among th em, we can 
mention the "Velocity Channel Analysis" (VCA, |Lazarian &| 
Pogosyan 2000; Esquivel et al. 2003, Lazarian & Pogosyan 



2004| |Lazarian et al.||2001[ |Padoan et al.||2003j |Chepurnov 
& Lazar ian 2009J), the "Spectral Correlation Function" (SCF 



Rosolowsky et al. 1999]|Padoan et aLl^OOlt 1 , "Velocity co- 
ordinate spectrum ( Lazarian & Pogosyan 2008 [ |2006), or 
"Modified Velocity Centro i ds" ([Lazarian & Esquivel |2003| 
Esquivel & Lazarian 2005; Ossenkopf et al. 2006; Esquivel 
et al.||2007[ |Chepurnov et al.||2006| |Ch epurnov & Lazarian 
2009| |Padoan et al.||2009| l Different methods are suitable in 
different situations, for instance VCA is adequate to retrieve 
velocity information from supersonic turbulence, while cen- 
troids are only useful in subsonic (or at most, transonic) tur- 
bulence. 

More sophisticated techniques for the analysis of the data, 
including using wavelets instead o f Fourier transforms ( |Hen-| 
riksen|[T994| |Stutzki et al.||1998} have been attempted! 
rather sophisticated example of such an alysis is a use of 
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the Pri ncipal Component A nalysis 2 (PCA) ( |Heyer & Schloerb 
[19971 |Brunt & Heyer|2002) . 

The techniques above provide the spectrum of turbulent ve- 
locity and/or density. Spectra, however, do not provide a com- 
plete d escription of turbulence. For example, Chepurno v et al.| 
( 2008 ) showed that fields with a very different distribution of 
density could have the same power spectrum. Intermittency 
and topology are some of the turbulence features that require 
additional statistical tools. Measures of turbulence intermit- 
tency include probability distributio n functions (PDFs, see for 
instance Falgaro ne et al.|1994]|Lis et al.|1996]|1998|[Klessen| 

1 The VCA and SCF uses the same sort of two point statistics. An insub- 
stantial difference is that VCA stresses the use of fluctuation spectra, while 
SCF uses structure functions. A more fundamental difference stem from the 
fact that the VCA is based on the analytical theory relating the observed spec- 
trum of intensity fluctuations and the underlying spectra of velocity and den- 
sity, while the choice of structure function normalization adopted in the SCF 
prevents the use of the developed theoretical machinery for a quantitative 
analysis I Lazarian 2009}. 

2 We should mention, however, that the use of integral transforms, other 
than the the Fourier transform, provide a new way of getting the same in- 
formation. For instance, using Fourier transforms one can also attempt to get 
the empirical relation between the underlying velocity and the statistics of the 
fluctuations. However, the VCA-obtained relations are very difficult to match 
empirically because of the strong non-linear nature of the transforms. There- 
fore it is not surprising that the empirically obtained relations for the PCA 
do not reproduce important cases discovered by the VCA. The PCA analysis 
failed to reveal the dependence of the Position-Position- Velocity (PPV) fluc- 
tuations on the density, which should be the case for the shallo w spectrum of 
fluctuations. See more on the testing the PPV dependencies in Che purnov &| 
Lazarian (20091. 
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2000 ), high er-order She-Leveque exponents (|She & Leveque 
JPadoan et al.||2003| |Falgarone et aLl|2005[ |Hily -Blant 
et al.| 2008). In~a ddition, genus (see also|Lazarian et al.|2002| 
Kim & Park 2007 ; Chepurno v et al.|2008| l was used to study 
the topological structure of the diffuse gas distribution. 

Astrophysical turbulence is very complex, e.g. it can hap- 
pen in the multi-phase media with numerous energy sources. 
Therefore the synergy of different techniques provide an in- 
valuable insight into the properties of the ISM turbulence 
and th e self-organization o f in terstellar medium. Fo r in- 
stance, |Kowal et aL] ( |2007| ) and |Burkhart efaL] ((2009) have 
applied several statistical measures, including a measure of 
bispectrum, to a wide set of MHD simulations and synthetic- 
observational data. More recent work applying the same set 
of techniques to the Small Magellanic Cloud (S MC) data re- 
vealed the advantages of such an approach (Burkhart et al. 
|2010| l. It is the consistency of the results obtained by differ- 
ent techniques, which made the study of turbulence properties 
in SMC more reliable. 

Thus it is important to proceed with th e quest for new tech- 
niques for astrophysical studies (see Lazarian 2009 and ref. 
therein). In this paper we add the Tsallis statistics to the ar- 
senal of tools available for ISM studies of turbulence. This 
is a description of the PDFs that have been successfully used 
previously in the context of solar wind observations and it 
shares similarities to previous studies of ISM turbulence (e.g. 
|Klessen|2000l >. 

The earlier solar wind relation between the Tsallis statistics 
and turbulence was sought on the basis of ID MHD simula- 
tions. We know that the properties of turbulence strongly de- 
pend on its dimensions. Thus we use 3D MHD simulations. 
We analyze a reduced set of MHD simulations, including col- 
umn density maps which are the most easily available data 
from observations. 

This paper is organized as follows. In section|2]we briefly 
present the Tsallis distribution, which we will use to char- 
acterize magnetohydrodynamic (MHD) simulations. The de- 
scription of the numerical models and the results (for 3D data 
as well as column density) can be found in sections [3j|5] We 
discuss the significance of our findings in|6] Finally, in section 
[7] we provide our summary. 

2. TSALLIS STATISTICS 



IBurlaga & Vifias| (|2004albl |2005a|b] ); Burlaga TNess &| 
Acuna ~i|2~006a); |Burlaga, Vinas, Ness, & Acuna| (|2006b ); 
Burlaga & Vinas (2006); Burlaga, Vinas, & Wang| p007b); 



Burlaga, Ness, & AcunaJ ^2007a| |2009|l , have studied the 
PDFs of the variation of the magnetic field strength in the 
solar wind, measured by the Voyager 1 and 2 spacecrafts. In 
particular, they analyzed the fluctuations of magnetic field in- 
crements [B(t + T m )-B(t)] as a function of the timescale r m , 
and showed that for a wide range of scales their PDFs could 



be de scribed by a Tsallis (or ^-Gaussian) distribution (Tsallis 
1988). The Tsallis PDF of increments (A/) of an arbitrary 
function has the form: 



R q =A 



„Af(t,r) 2 
l + (q-l) — 



-1/(9-1) 



(1) 



where q is the so called "entropic index" or "non-extensivity 
parameter" and is related to the size of the tail of the dis- 
tribution, w provides a measure of the with of the PDF (re- 
lated to the dispersion of the distribution), while A measures 



the amplitude. Notice that the Tsallis distribution is symmet- 
ric 3 . Varying the parameter q in the Tsallis distribution one 
obtain distributions that range from Gaussian to "peaky" dis- 
tributions with large tails. The parameter q is closely related 
to the kurtosis (fourth order one-point moment) of the PDF, 
and similarly the parameter w is related to the variance of the 
PDF. However, they are not the same, while fitting a Tsallis 
distribution to a PDF, the two parameters are calculated si- 
multaneously, therefore w will depend on the value of q and 
vice versa, while the dispersion of a distribution is indepen- 
dent of the kurtosis. In the next section we discuss how this 
difference might play in our favor. Furthermore, the Tsallis 
distribution is often used in the context of n on-extensiv e sta- 
tistical dynamics, it was originally derived (Tsallis 1988} from 
an entropy generalization to extend the traditional Boltzmann- 
Gibbs statistics to multifractal systems (such as the ISM). The 
Tsallis distribution reduces to the classical Boltzmann-Gibbs 
(Gaussian) distribution in the limit of q — ► 1. It provides a 
physical foundation to other functions that have been used 
to model empirically velocity PD Fs in space plasmas (such 
as the kappa function, see Maksim ovic et al.||1997[ [Leubner 
' ' However, the physical interpretation of ISM simula- 



tions/observations PDFs in the light of such statistical dynam- 
ics is beyond the scope of this paper. 

Turbulence intermittency manifests itself in non-Gaussian 
distributions, thus one can hope to establish a link between 
the properties of turbulence and the Tsallis distributions pa- 
rameters obtained from simulations and/or observations. This 
provides an alternative/complementary method to higher or- 
der statistics (Burkh art et al.|2009|l, or to other PD F analysis 
(such as those in Falgaro rie et al.| l994|JKlessen|2000 where 
the kurtosis and/or skewness of the distributions is measured), 
all of which are sensitive to the non-Gaussian nature of real 
turbulence. 

Following the procedure presented by Burlaga and collabo- 
rators in what follows we make histograms of the distribution 
of fluctuations from our simulations. But, in contrast to their 
studies, where only the magnetic field variations were studied 
at different times, we consider the spatial variations of density, 
velocity and magnetic field intensity at a given time, as a func- 
tion of a spatial scale (or lag) r. To make the fits and present 
the results we remove the mean value of the increments and 
scale them to units of the standard deviation of each PDF. In 
other words, we have used A/(r) = [f(r)— (/(r)) x ]/cr/, with 
f(x) = p(x+r)-p(x), v x (x+r)-v x (x), v Y (x+r)-Vy(x), v z (x+r)- 
v z (x), B x (x + r)-B x (x), B y (x + r)-B y (x), and B z (x + r)-B : {x), 
where (...) x stands for spatial average. We have n eglected 
the anisotropy introduced by the magnetic field (see |Cho &| 
Lazarian 2009 ), the lag r is considered along each of the three 
cardinal directions and the resulting histograms include vari- 
ations in all of them. 

3. MHD SIMULATIONS 

We have taken a reduced subset of th e MHD s imulations 
of fully-developed (driven) turbulence of |Kowal et al.| (|2007). 
We use the data output of the simulations to study the result- 
ing PDFs of density, magnetic field and velocity in different 
turbulence scenarios. 

The code used (described in detail in Ko wal et al.| |2007) 

3 Measurements of velocity increments in t he solar wind show asymmetric 
PDFs, which has lead Burlaga & Vinas 1 2004a) to add a cubic term to the ex- 
pression in eq. l[TJ, in order to lorm a generalized" (non-symmetric) Tsallis 
distribution that can be used to describe skewed PDFs. 
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TABLE 1 
Parameters of the simulations 



Model 


P 




AV 


M A b 


Description 


1 


0.1 


0.1 


2.0 


2.0 


Supersonic & super-Alfvenic 


2 


0.1 


1.0 


2.0 


0.7 


Supersonic & sub-Alfvenic 


3 


1.0 


0.1 


0.7 


2.0 


Subsonic & super-Alfvenic 


4 


1.0 


1.0 


0.7 


0.7 


Subsonic & sub-Alfvenic 



a Sonic Mach number, defined as Ms = (v/cs) (averaged over the entire 
computational domain). 

b Alfven Mach number, defined as Ma = ( v /va) (averaged over the en- 
tire computational domain). 

solves the ideal MHD equations in conservative form: 

dp 



dt 



+ V-( j ov) = 0, 



IF 



+v- 



pv\ + p + 



<9B 

Ik 



I 



1 



BB 



= f. 



47T 

Vx(vxB) = 0, 



(2) 
(3) 
(4) 



along with the additional constraint V • B = in a Cartesian, 
periodic domain. An isothermal equation of state p= c 2 s p is 
used, where p is the mass density, v the velocity, B the mag- 
netic field, p the gas pressure, and c s the isothermal sound 
speed. The term f in the equation of momentum (3} is a 
large scale driving. The driving is purely solenoidal, ap- 
plied in Fourier space at a fixed wave-number k = 2.5 (a scale 
2.5 times smaller than the size of the computational box). 
The rms velocity is maintained close to unity, so that v is in 
units of the rms velocity of the system, the Alfven velocity 
va = B/(47rp) 1 / 2 is also in the same units. The magnetic field 
is of the form B = B ex t + b, that is a uniform background field 
B ext plus a fluctuating part b (initially equal to zero). As it is 
customary, the units are normalized such that the Alfen veloc- 
ity va = B e xt/(47rp) 1 / 2 = 1 and the mean density po = 1. This 
way, the controlling parameters are the sound, and the Alfven 
speeds (i.e. the gas pressure and the magnetic pressure, re- 
spectively), yielding any combination of subsonic or super- 
sonic, with sub-Alfvenic or super-Alfvenic regimes of turbu- 
lence. For this paper we explore the different combinations of 
parameters summarized in Table [T] 

We should note that it has been recently found that a com- 
pressive component of the driving can have an important ef- 
fect on the statistical properties of the flow (|Schmidt et al.| 
[2008) [Feden-ath ~eTal.|2008[ |Schmidt et al.|2009| l. This is an 
interesting effect that require further study, but for the time 
being we will restrict ourselves to the usual divergence-free 
driving. 

4. STATISTICS OF THREE-DIMENSIONAL DATA 

Below, we analyze the results of our simulations to test 
whether the statistics of the PDFs of the fluctuations of the 
variables can be described by the Tsallis formalism. To to 
this we fit to each histogram a Tsallis distribution with the 
Levenberg-Marquardt algorithm ( |Press et al. 1992), the his- 
tograms and the fits are presented in Figure [l] the symbols 
are the data from the simulations, and the lines are the corre- 
sponding fits. For visual purposes we have shifted vertically 
the results by successive factors of 100 for each diffe rent lag- 
Simila rly to the results from solar wind o bservations (Burlaga 
& Vinas 2004a b 2005a b ; Burlaga et al. 2006a b ; Burlaga & 

I _l I 1 1 I 1 1 2 I I 1 1 S 



Vifias|2006| [Burlaga et d.|2007a|b||2009) , the magnetic field 
components, as well as the density and velocity components, 
show kurtotic PDFs that become smoother as we increase the 
scale-length (lag). One difference is the fact that the velocity 
increments in our simulations are symmetric, as opposed to 
the skewed distributions observed in the solar wind at small 
scales. It is also interesting that (only) for supersonic turbu- 
lence density distributions are quite distinct from other quan- 
tities. Such behavior resembles the power spectrum, where 
the magnetic field and velocity for supersonic turbulence have 
steep spectra regardless of sonic Mach number, while density 
becom es shallower as the Mach nu mber increases (see for in- 
stance Esquivel & Lazarian 2005). Such shallower slopes of 
the density power-spectrum have been also reported when the 
driving changes from solenoidal to c ompressive (potential ) 
for a fixed (> 1) sonic Mach number (Federrath et al. 2009). 
This means that shocks in supersonic turbulence, which create 
more small-scale structure in density, yield not only shallower 
spectra, but also a larger degree of intermittency. At the same 
time, such shocks can be formed easier if the turbulence driv- 
ing mechanism has compressive modes. 

We can observe from Figure [T]that the Tsallis distributions 
fit remarkably well the PDFs in the simulations. Only small 
departures are evident at the tails of the PDFs, in particular 
for density in models 1 and 2 (supersonic turbulence). The 
values of the q and w parameters from the fits are plotted as a 
function of the separation r in Figures [2j and [3] respectively. 

In Figures [2] and [3] we see the same general trends found 
in solar wind st udies (Burlag a & Vifia s 2004a b~| |2005a|b| 
|Burlaga et al.|2006a|bj|Burlaga & Vinas|2006| |Burlaga et aL 



2007a b 2009) and electron MHD simulations (Cho & Lazar- 
|ian|2009| : a q parameter that decreases, along with a w value 
that increases as we take larger separations. In other words, 
the PDFs become more kurtotic, and narrower, as we go to 
small scales, while at large scales the PDFs are more Gaus- 
sian and wider. The same behavior (i.e. a decreasing kurtosis 
for an increasing lag) has been found also for velocity cen- 
troids in simulations and observations (see Falgarone et al. 



[19941 |Lis et al.|1998||Klessen|2000l|Federrath et al.|20 09). In 
fact, we have computed the kurtosis and variance of the PDFs 
(not shown here) and found a similar monotonic decrease and 
increase, respectively. Both the kurtosis and q tend to their 
Gaussian value (3 and 1, respectively) as the separation tends 
to the injection scale. 

It is noticeable that the density PDFs in the supersonic 
regime have a very distinct shape than those of the velocity 
and magnetic field components. In the models of supersonic 
turbulence they are much narrower, and more peaked than in 
subsonic turbulence. They also show a shallower slope in the 
q and w vs r plots. We can also notice slightly larger values 
of q (and lower w) for the velocity components with respect 
to the magnetic field components, more pronounced for high 
Alfen Mach numbers. 

5. STATISTICS OF OBSERVABLE DATA: COLUMN 
DENSITY 

The measures of non-Gaussianity presented in the previous 
section are interesting in their own right, and might provide 
us with insights to the intermittency of MHD turbulence from 
simulations. However, observations of the ISM cannot pro- 
vide direct three-dimensional information of the density, ve- 
locity, and magnetic fields. For instance, from spectroscopic 
observations we can access only the velocity distribution of 
the emitting material at a given position (or positions) in the 



4 



Esquivel & Lazarian 



(a) M s =2.0, M A =2.0 (b) M s = 2.0, M A =0.7 (c) M s = 0.7, M A =2.0 (d) M s = 0.7, M A =0.7 




f(x + r)-f(x) 



FIG. 1 . — PDFs of fluctuations for different spatial lags r, for all the models. The symbols represent data from the numerical simulations and the lines are the 
fits with a q-Gaussian (Tsallis) distribution. The different MHD variables are color coded according to the legend at the upper left corner, and the value of the lag 
is indicated at the left of each set of curves. 



plane of the sky. 

With the simulations presented in the previous section we 
have constructed (2D) maps of column density, assuming that 
the emitting material is optically thin, and with an emissiv- 
ity linearly proportional to the density (e.g. H I). From the 
two-dimensional column density maps we constructed PDFs 
of column density increments, and following the same pro- 
cedure we described before, we made fits to Tsallis distribu- 
tions. The PDFs and the fits are presented in Figure H] 

From the figure we can confirm the same general behavior 
of the PDFs, which become more Gaussian as we take larger 
separations. We can see that the PDFs are still consistent with 
Tsallis distributions, although the dependence on the lag is 
more subtle than with the 3D fields used in the previous sec- 
tion. 

In Figure [5] we present the parameters q and w as obtained 
from the fitting procedure. Both panels of Figure [5] confirm 
the decrease of q along with an increase of the width of the 
PDFs (yv) with the lag. From the figure it is also evident 
that the fitting parameters can clearly distinguish between 
the supersonic and subsonic models. However, this can be 
done with other mea sures, for instance the variance of density 
or column density (|Padoan et al.|[l997 |Passot & Vazquez- ] 
|Semadeni|2003l|Federrath et al.|2008||2009) , or t he skewness 
and kurtosis of density and column density PDFs (Kowal et al. 



|2007||Burkhart et al.|2009] ). A more difficult task is to extract 
information of the Alfven Mach number. One can see from 
Figure B] that the Tsallis parameters depend strongly on the 
sonic Mach number, and less so on the Alfven Mach num- 
ber. However, an encouraging clear dependence on the Alfven 
Mach number is evident in the supersonic models, in partic- 



ular for w (lower panel of Figure [5]), where all the Ai A = 2 
points lie above the M.a = 0.7 ones. But, for the subsonic 
models (M. s = 0.7) it not clear at all, and the behavior with 
the Alfen Mach number is the opposite, higher w values for 
sub-alfvenic models. 

6. DISCUSSION 

We found that the Tsallis expression [see Eq. ([T}] describes 
well the statistics of the PDFs of our 3D-MHD numerical sim- 
ulations. This is welcome news for both the ISM, where we 
propose to use Eq. ([T]i as a new tool to study turbulence, and 
for the solar wind studies, where the Tsallis expression was 
used, but its numerical testing was limited. 

We have to add that the present work is intended only as a 
proof of concept. While the Tsallis PDFs do provide physical 
grounds to interpret the distributions obtained in MHD simu- 
lations and/or observations, it is not clear whether the formal- 
ism of non-extensive statistical dynamics is applicable to the 
turbulent ISM. This, we feel deserve a separate study. How- 
ever, the fact that the expression provides both the fit to nu- 
merical MHD simulations and experimental data, and more- 
over seems to distinguish the Alfven Mach number in some 
cases, is encouraging. 

In regard to observational quantities, we should note that 
PDFs of velocity centroids (as opposed to column density, 
which is what we ave used in the present paper) have been 
employed t o characterize the inter mittency both i n numeri- 
cal models ( Falgarone et al.|[l994 Klessen 2000 Federrath 
|et al.|2009] > and molecular c loud observations (for instance p- 
Ophiu chi in Lis et al.|1998| Polaris and Taurus in |Hily-Blant 
|et al.|2008) . There is always a trade off while choosing dif- 
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(a) M s = 2.0, M A =2.0 



(b) M s =2.0, M A =0.7 




FIG. 2. — Values of the q parameter from the fits shown in Figure[T| the symbols (and colors) indicate the MHD variable used, according to the legend at the 
left of the plots. 



ferent quantities to measure turbulence. On the one hand ve- 
locity is a dynamical quantity that is easy to relate with turbu- 
lence theories, however (at least in the case of power-spectra) 
is trickier to obtain, and in fact (for power-spectra) velocity 
centroids fail to trace the velocity statistics in strongly super- 
sonic turbulence (see for instance Esquivel & Lazarian 2005 ). 
On the other hand density fluctuations are not only produced 
from turbul ence, but are more robustly o btained from obser- 
vations (e.g. Lazarian & Pogosyan 2008). We must stress that 
the healthy approach is to use these tools, and others as well, 
as complementary. 

The most important finding of this paper is the dependence 
of the q and w parameters entering Eq. (fT) on both sonic and 
Alfven Mach number of the turbulence!~These numbers are 
essential f or many astrophysical pr ocesses from star forma- 
tion (see McKee & Ostriker 2007), magnetic diffusion en- 
abled by reconnection (Lazarian 2005), propagation of heat 
(Nara yan & Me dvedev 2 001 [ |Laza rian 2006a}, co smic rays 
( Yan & Lazarian 2004 ) and dust grains acceleration dYan et al. 
2004). 

By now several ways of estimating the Alfven and sonic 
Mach numbers have been proposed (see Kowal et al. 2007), 
all with their own limitations and uncertainties. Thus the new 
approach here based on Eq. ([T]) is a welcome addition to the 
existing tools. 

We should add however, that a number of turbulence pa- 
rameters, which are as important as these Mach numbers, are 
not being explored in this work. The scale of turbulence, for 
instance, has important consequences for sta r formation (e.g. 
|Klessen et al.|2000[|Mac"Low & Klessen|2004) as well as for 



dynamo theory (e. g. Brandenburg & Subramanian 2 005b|a 
|Hanasz et al.|2009) 

7. SUMMARY 

We have used three-dimensional MHD simulations to study 
the distribution functions of increments of density, velocity, 
and magnetic field. In addition, we studied synthetic column 
density maps, which can be directly compared with observa- 
tions. The results can be summarized as follows. 

1 . The PDFs of each of all the quantities is well described 
by a Tsallis distribution. Such distributions are sym- 
metric, and characterized basically by three parameters: 
an amplitude A, a width w , and the entropic index q that 
relates to the shape (large values yield peaked distri- 
butions, while with small values they approximate to a 
Gaussian). 

2. The general trend of the PDFs for different increments 
are consistent with previous studies, of solar wind ob- 
servations. That is, the PDFs are close to a Gaussian 
distribution for large separations, and become more 
kurtotic and narrower (i.e. a larger degree of intermit- 
tency) for small separations. 

3. While the PDFs of (3D) density are quite similar to 
those of magnetic field and velocity in subsonic turbu- 
lence (regardless of Alfven Mach number), for super- 
sonic turbulence (also independently of A4a) they are 
noticeable narrower and more peaked. That is, the fits 
of the PDFs to a Tsallis distribution have smaller w, and 
larger q parameters. 
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FIG. 3. — Values of the w parameter from the fits shown in Figure[T] the symbols (and colors) indicate the MHD variable used, according to the legend at the 
left of the plots. 



The column density maps show smoother distributions 
compared to those of density in 3D, however, we found 
that the q values are larger for supersonic turbulence 
(and lower w as well), regardless of the Alfven Mach 
number. 

The w parameter of the Tsallis fits to the PDFs of col- 
umn density increments seems to show a sensitivity to 
the Alfven Mach number, at least in the cases of super- 
sonic turbulence. With the limited data set used here is 
hard to conclude if this behavior is general, and whether 
it is useful to study the magnetization of ISM turbu- 
lence. Clearly, future studies of the PDFs in ISM turbu- 
lence (with a wider parameter coverage) are necessary. 
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FIG. 5. — Parameters q and w from d the Tsallis distribution fits shown in 
Figure|4] 



